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It is generally accepted that statistics of energy levels in closed chaotic quantum systems is 
adequately described by the theory of Random Hermitian Matrices. Much less is known about 
properties of "resonances" - generic features of open quantum systems pertinent for understanding 
of scattering experiments. 

In the framework of the Heidelberg approach to quantum chaotic scattering open systems are 

characterized by an effective non-Hermitian random matrix Hamiltonian H^j . Complex eigenvalues 

of _ffe/ are S-matrix poles (resonances). We show how to find the mean density of these poles 

t~^ ' (Fyodorov and Sommers) and how to use the effective Hamiltonian to calculate autocorrelations of 

0^ \ the photodissociation cross section (Fyodorov and Alhassid). 

CTS . In the second part of the paper we review recent results (Fyodorov, Khoruzhenko and Sommers) 

on non-Hermitian matrices with independent entries in the regime of weak Non-Hermiticity. This 
^ , regime describes a crossover from Hermitian matrices characterized by Wigner-Dyson statistics of 

D ■ real eigenvalues to strongly non-Hermitian ones whose complex eigenvalues were studied by Ginibre. 
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I. INTRODUCTION 



^ ' As is well-known, statistics of highly excited bound states of closed quantum chaotic systems of quite different 



microscopic nature is universal. Namely, it turns out to be independent of the microscopic details when sampled on 

the energy intervals large in comparison with the mean level separation, but smaller than the energy scale related by the 

i^ , Heisenberg uncertainty principle to the relaxation time necessary for the classically chaotic system to reach equilibrium 

,_^ in the phase space [|l|. Moreover, the spectral correlation functions turn out to be exactly those which are provided by 

f — . _ the theory of large random matrices on the local scale determined by the typical separation A = {Xi — Xi^i) between 

On ■ neighboring eigenvalues situated around a point X , with brackets standing for the statistical averaging 0. Microscopic 

^ ' justifications of the use of random matrices for describing the universal properties of quantum chaotic systems have 

^^. been provided recently by several groups, based both on traditional semiclassical periodic orbit expansions ||3|,y and 

'^ ' on advanced field-theoretical methods [pp| . These facts make the theory of random Hermitian matrices a powerful 

J^ , and versatile tool of research in different branches of modern theoretical physics, see e.g. ||,|j. 

ii ' Very recently complex eigenvalues of non-Hermitian random matrices have also attracted much research interest 

r^ ] due to their relevance to several branches of theoretical physics. Most obvious motivation comes from the quantum 
O ' description of open systems ||8|^|lC| whose fragments can escape, at a given energy, to infinity or come from infinity. 
^ ■ For systems of this kind the notion of discrete energy levels loses its validity. Actually, chaotic scattering manifests 
itself in terms of a high density of poles of the scattering matrix placed irregularly in the complex energy plane. Each 
of these poles, or resonances, Ek — £k — fTfc, is characterized not only by energy £k but also by a finite width Ffe 
defined as the imaginary part of the corresponding complex energy and reflecting the finite lifetime of the states in 
the open system. Recently, the progress in numerical techniques and computational facilities made available high 
accuracy patterns of resonance poles for realistic atomic and molecular systems in the regime of quantum chaos, see 
e.g. |l|]I|. 

Due to the presence of these resonances, elements of the scattering matrix show irregular fluctuations as functions 
of the energy of incoming waves, see |Tq] and references therein. The main goal of the theory of quantum chaotic 
scattering is to provide an adequate statistical description of such a behavior. 

Whereas the issue of energy level statistics in closed chaotic systems was addressed in an enormous amount of 
papers statistical characteristics of resonances are much less studied and attracted significant attention only recently, 
see [p|-[ll| and references in |16| . 
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One possible way of doing this is to address resonances in the so-called " Heidelberg approach" suggested in the 
pioneering paper |171| and described in much detail in |16| . The approach turns out to be the most natural framework 
for incorporating random matrix description of the chaotic scattering. 

The starting point of this approach is a division of the Hilbert space of the scattering system into two parts: 
the "interaction region" and the "channel region". The channel region is supposed to describe a situation of two 
fragments being apart far enough to neglect any interaction between them. Under these conditions their motion 
along the collision coordinate is described by a superposition of incoming and outgoing plane waves with wavevectors 
depending on the internal quantum states of the fragments. We assume that at given energy E exactly M different 
quantum states of the fragments are allowed, defining M "scattering channels" numbered by the index a. 

At the same time, the second part of the Hilbert space is to describe the situation when fragments are close to 
one another and interact strongly. Correspondingly, any wavefunction of the system \^{E)) can be represented as 

two-component vector: |$(£')) = , J, with u and ip describing the components of the wave function inside the 

interaction (respectively, channel) region. 

Using the standard methods of the scattering theory exposed in detail in llq] one can relate two parts of the 
wavefunction to one another and finally arrive at the following representation of the energy-dependent scattering 
matrix S in terms of an effective non-Hermitian Hamiltonian TCef = H ~ iT: 

SabiE) = Sab - 2tTr Y,Wa^[E - Hefj^'W.b (1) 

ij 

with the Hermitian Hamiltonian H describing the closed counterpart of the open system (i.e. interaction region 
decoupled from the channel one) and the anti-Hcrmitian part T arising due to a coupling to open scattering channels. 
In this expression the Hamiltonian H is written in some arbitrary basis of states |i), such that Hij = {i\H\j). The 
amplitudes Wai, a, — 1,2, ...,M are matrix elements coupling the internal motion in an "internal" state \i) to one 
out of M open channels a. One also has to choose the anti-Hermitian part to be F = '"'X^a ^iaW*^ in order to ensure 
the unitarity of the M x M scattering matrix S{E) ||l^-po|. 

It is natural to expect, that universal properties of open chaotic systems are inherited from the corresponding 
universality of levels of their closed counterparts. Of course, one can expect a relation of this kind only when 
incoming particles stay inside the interaction region long enough to be able to experience the chaoticity of internal 
dynamics. Going from the time domain to the energy domain, this fact suggests that only scattering characteristics 
on a scale shorter than inverse classical relaxation time on the energy shell are expected to be universal. Another 
characteristic energy scale in this domain is a typical level spacing A of the closed counterpart of our quantum open 
system. Thus, we expect the scattering characteristics (in particular, the statistics of resonances) to be universal on 
the scale comparable with A. In contrast, smooth energy dependence of S"— matrix elements on a much larger energy 
scale must be system-specific. 

The next step is to incorporate the random matrix description of quantum chaotic systems by replacing the Hamilto- 
nian H hy a, random matrix of appropriate symmetry. Namely, chaotic systems with preserved time-reversal invariance 
(TRI) should be described by matrices Hij which are real symmetric. Such matrices form the Gaussian Orthogonal 
Ensemble, whereas for systems with broken TRI one uses complex Hermitian matrices from the Gaussian Unitary 
Ensemble [|§. 

The third essential ingredient of the Heidelberg approach is performing the ensemble averaging non-perturbatively 
in the framework of the so-called supersymmetry technique. It was invented initially by Efetov in the context of 
theory of disordered metals and the Anderson localization [pl| ,p2[ and adjusted for the description of quantum chaotic 
scattering by Verbaarschot, Weidcnmiiller and Zirnbauer |17| . 

The Heidelberg approach turns out to be a very powerful tool for extracting universal properties of open chaotic 
systems. In the first part of the paper we outline the derivation of the mean resonance density in the complex plane 
following Fyodorov and Sommers P| . Another recent example of the utility of the non-Hermitian effective Hamiltonian 
H — iT is that its resolvent is related to the probability for an excited system to decay via open channels. For this 
reason it can be used to calculate the statistics of such quantities as e.g. photodissociation cross-section in the regime 
of quantum chaos. We present the formula for cross-section autocorrelation function derived recently by Fyodorov 
and Alhassid (2|]. 

The fact that non-self adjoint operators appear quite generally when one considers open systems of various types 
is known for a long time pO|. It is therefore not surprising that open quantum systems were the first examples of 
applications of non-Hermitian random matrices |p|, 3Hl6|, see also recent papers p0| , p^ , p5[ . Other early applications 
included also studies of dissipative quantum maps |26H28|] and chaotic dynamics of asymmetric neural networks [E9|. 



Recently, however, random matrices (more generally, random linear operators) with complex eigenvalues emerged 
in many other physical contexts. Let us briefly mention the most interesting examples. 

• The effective Hamiltonian describing a thermal motion of an isolated vortex in disordered type-II supercon- 
ductors with columnar defects has a form of that for a quantum particle in an imaginary vector potential 
A: 

H = ^{~ih\/ + iAf + V{r), (2) 

/TO 

with V{r) being a random potential generated by defects. The imaginary vector potential makes the Hamiltonian 
to be a non-Hermitian one. This fact pointed out by Hatano and Nelson [Q gave a boost to several interesting 
studies |l|,|2|. 

• A classical diffusing particle advected by a stationary random velocity field v is described by a non-Hermitian 
Fokker-Plank random operator Cpp @| : 

^n(r, t) - CFp72ir, t) = {DV^ - Vv) n(r, t) (3) 

where n(r,i) is the density of particles and D is the diffusion constant. 

• Recent attempts to understand the universal features of chiral symmetry breaking in Quantum Chromodynamics 
required to consider quarks in a finite chemical potential /i interacting with the Yang-Mills gauge field. The 
corresponding partition function is given by: 

Z{m,^l) = ^'Y{Aei(mJ+^i%+'Dj\ (4) 

where V = J^d^ + i^fiAfi is the Euclidean Dirac operator in the gauge vector potential A, and (. . .)_^ = 
J DA{. . .)exp— J d'^xF^j^ stands for the averaging over all configurations of the gauge field F^^ = d^A,y — 

A finite chemical potential fi makes the corresponding operator to be a non-selfadjoint one with complex eigen- 
values. This fact makes a problem of numerical evaluation of the partition function by lattice simulations to 
be a very difficult one. Recently it was suggested, that some universal features of the model can be correctly 
recovered if one replaces the true gauge-field averaging by averaging over random gauge- field configurations. As 
a result, one comes to a class of non-Hermitian random matrix problems of a particular type [ p4| . 

• Finally, let us mention that there exist several interesting links between complex eigenvalues of non-Hermitian 
random matrices and systems of interacting particles in one and two spatial dimensions ||35| , |36| . 

At the same time, our knowledge of statistical properties of random non-selfadjoint matrices is quite scarce and 
incomplete. This fact recently stimulated efforts of different groups to improve our understanding in this direction 

0- il. 

Traditional mathematical treatment of random matrices with no symmetry conditions imposed goes back to the 
pioneering work by Ginibre Jl3[ who determined all the correlation functions of the eigenvalues in an ensemble of 
complex matrices with independent Gaussian entries. The progress in the field was rather slow but steady, see 
~|,||J29[g|||,|^,|9|. 



Surprisingly, all these studies completely disregarded the existence of a nontrivial regime of weak non-Hermiticity 
recognized in the work by Fyodorov, Khoruzhenko and Sommers Bq] , see more detailed discussion in [ pO[ . The guiding 
idea to realize the existence of such a regime comes from the experience with resonances ^ . Guided by that example 
one guesses that a new regime occurs when the imaginary part of typical eigenvalues is comparable with a mean 
separation between neighboring eigenvalues along the real axis. 

One can again employ the same supersymmetry approach which was used to study resonances and obtain the mean 
density of complex eigenvalues in the regime of weak non-Hermiticity for matrices with independent elements |^,Q . 
The density turned out to be described by a formula containing as two opposite limiting cases both the Wigner 
semicircular density of real eigenvalues typical for Hermitian random matrices and the uniform density of complex 
eigenvalues discovered for strongly non-Hermitian random matrices already by Ginibre B3] , in much details addressed 
by Girko Ea] and studied for different cases by other authors p3,Ea,E9| . 



Very recently, Efetov p2| showed the relevance of almost-Hermitian random matrices to the very interesting problem 
of motion of flux lines in superconductors with columnar defects |30|] . He also managed to derive the density of complex 
eigenvalues for a related, but different set of almost-symmetric real random matrices. This development clearly shows 
that, apart from being a rich and largely unexplored mathematical object, almost-Hermitian random matrices enjoy 
direct physical applications and deserve to be studied in more detail. 

Actually, the non-Hermitian matrices considered in |3q ] and [ ||32| ] are just two limiting cases of a general three- 
parameter family of non-Hermitian ensembles [ pO| . In second section of the paper we outline the derivation of this 
fact and present the resulting expression in terms of a non-linear cr-model integral. 

Although giving an important insight into the problem, the supersymmetry non-linear a— model technique suffers 
from at least two deficiencies. The most essential one is that the present state of art in the application of the 
supersymmetry technique gives little hope of access to quantities describing correlations between different eigenvalues 
in the complex plane due to insurmountable technical difficulties. At the same time, conventional theory of random 
Hermitian matrices suggests that these universal correlations are the most interesting features. The second drawback 
is less important for a physicist, but a crucial one for the mathematicians: at the moment the supersymmetry technique 
can not be considered as a rigorous mathematical tool and has the status of a powerful heuristic method. 

Fortunately, for the simplest case of almost-Hermitian Gaussian random matrices one can develop the rigorous 
mathematical theory based on the method of orthogonal polynomials. Such a method is free from the above mentioned 
problem and allows one to study correlation properties of complex spectra to the same degree as is typical for earlier 
studied classes of random matrices |3^. We briefly discuss the obtained results in the end of the paper. The detailed 
exposition of the method and the derivation of the results can be found in ||5^ . Unfortunately, the paper Q contained 
a number of misleading misprints. For this reason we indicate those misprints in the present text by using footnotes. 

II. NON-HERMITIAN RANDOM MATRICES IN THE THEORY OF CHAOTIC QUANTUM 

SCATTERING 

To calculate the density of resonance poles in the complex energy plane we notice that they are just eigenvalues of 
the effective non-Hermitian Hamiltonian introduced in Eq.(nf). 

Without any loss of generality coupling amplitudes Waj can be chosen in a way ensuring that the average 5— matrix 
is diagonal in the channel basis: Sab = SabSaa- Then one finds the following expression pTJ ]: 

^--^ l+7,g(X) ' "f-=''2^^^aW^a (5) 



where g{X) = iX/2 + 7ri^(X) and tti>{X) = (1 — Ar^/4)^/^ is the semicircular level density. The strength of coupling 
to continua is convenient to characterize via the transmission coefficients Ta — 1 — \Saa\'^- These quantities measure a 
part of the incoming flux in a given channel that spends a substantial part of the time in the interaction region ||l^,plj . 
The case T^ <C 1 corresponds to almost closed channel a, whereas the opposite limiting case Ta = 1 corresponds to 
the perfect coupling between the interaction region and the channel a. It is easy to see that both limits 7a — > and 
7a — > c» equally correspond to the weak effective coupling regime whereas the strongest coupling (at fixed energy X 
) corresponds to 7a = 1. 

In the case of weak effective coupling to continua individual resonances do not overlap: (F) <C A, with A = 
{i'{X)N)^^ standing for the mean level spacing of the "closed" system and (F) standing for the mean resonance 
width. Under these conditions one can use a simple first order perturbation theory to calculate resonance widths in 
terms of eigenf unctions of the closed system. Since different components of eigenvectors of large random matrices are 
decorrelated and Gaussian-distributed, one finds in such a procedure that the scaled widths ys = tSv are distributed 

according to the so-called x^-distribution: 

where F(z) stands for the Gamma function and the parameter p = M {v = 2M) for systems with preserved (broken) 
time reversal invariance, M being the number of open scattering channels. The case z^ = 1 is known as Porter-Thomas 
distribution and was shown to be in agreement with experimental data (see some references in |52|]). 

Experimentally, one quite frequently encounter the case of only M ^ 1 open channels and (F) ~ A, see e.g. [p3[-p5[. 
Under this situation resonances overlap considerably and one can not use perturbation theory any longer. The 
problem of determining the statistical characteristics of the chaotic scattering in the regime of (partly or completely) 



overlapping resonances is of essentially non-perturbative nature. As a result, one has to use some non-perturbative 
methods allowing to evaluate the ensemble averaging exactly for arbitrary ratio (r)/A. 

Fortunately, one can study very efficiently various universal statistical features of chaotic quantum scattering by 
performing the ensemble averaging with the use of the supersymmetry method |17| . 

One can recover the spectral density 

N N 

p{Z) = J2 ^^'^ {Z - Zu) = Y.5{X ^ X,)S{Y - Y,) EE p{X, Y) (7) 

k=l fc=l 

of complex eigenvalues Z^ = X^ + iYk, A: = l,2,...,iVif one knows the "potential" |E9[: 



$(X,y, n)^^ \nDet[{Z - nef){Z - HefV + k^] 
Ztt 

in view of the relation:p(X, F) = liiaii^^od^ ^{X,Y,k), where d^ stands for the two-dimensional Laplacian. Techni- 
cally, it is convenient to introduce the generating function (cf. |16|] ) 

^ Det[{z-nef){z-nefy + K^] 

Det [{Zb - TiefKZb - He/)t + k2] ^ > 

in terms of which 

p(Z) = --lim-^ lim -^Z (9) 

^^ ' TT «^o dZ* z,^z dZb 

To facilitate the ensemble averaging we follow the standard route and represent the ratio of the two determinants 
in Eq.(p|) in terms of a Gaussian superintegral over eight-component supervectors $i = I j^, I ] where ^i(±) = 

^'(±) ^ and 

with components ri{+),ri(—): i ~ 1,2,..., N being complex commuting variables and Xi(+)iXi(^) forming the 
corresponding Grassmann parts of the supervectors ^^ (±). 

Further evaluation goes along the lines discussed in |1^ in more detail. After a set of standard manipulations 
one arrives at the following expression for the density px{y) = ^p{X:Y)A'^{X) of scaled resonance widths y = ^ 
(measured in units of the local mean level spacing A(X) of the closed system) for the resonances whose positions are 
within a narrow window around the point X of the spectrum: 

. M r . - 

{px{y)) - Y^ y d^i{Q)Stv (af )q) Str (a,g) exp ^yStr (a.g) J] Sdet^^/" 1 - ^ {Q,a,} (10) 

Here the integration goes over the set of 8 x 8 supermatrices Q satisfying the constraint Q^ — ~1 and ] QjO't f = 

Q&r + a-rQ stands for the anticommutator. Properties of these matrices and the integration measure dfi{Q) can be 
found in [2^. Other 8x8 supermatrices entering the expression Eq.(|o|) are as follows: 

fa 
f3 

and fa; , f 3 ' are 4x4 diagonal supermatrices: fa = diag{f, f}; fg ' = diag{02,f}. with f = diag(l,— 1). We also 
introduced quantities ga = 2;;^nn('>'a + 70"^) related to the transmission coefficients as ga = 2/Ta — 1 and used the 
symbols Str, Sdet for the graded trace and the graded determinant, correspondingly. 

The expression above is valid for chaotic systems with preserved as well as with broken time-reversal invariance. To 
extract the explicit form of the distribution function one still has to perform the integration over the manifold of the 
supermatrices Q which is different for two cases. In general it is a rather difficult calculation due to a cumbersome 




parametrisation of that manifold. At the moment the result is known for the simplest case of the systems with 
broken time-reversal invariance PJigI . For the sake of simplicity we present this distribution for the case of equivalent 
channels a = 1, ..., M with equal transmission coefficients Ta = T . 

First of all, it turns out that the mean resonance width is related to the transmission coefficient T as: 

(r>^-A|lln(l-T).A|^ln^ (11) 

The formula Eq.(ni]) is well known in nuclear physics as Moldauer-Simonius relation | |5q |. 

It is convenient to use the parameter k = — -^ In (1 — T) as a measure of the resonance overlap. It is related to the 

mean widths as k = tt-M- and therefore gives a typical number of neighboring resonances that overlap substantially. 
Measuring the resonance widths in units of the mean widths (F) one finds the following distribution function. 

P ivs - 7^) = .r^^L . r^ dtt^'^-' (12) 

'Vsh- 



(F)y 2F(Af)K2/2 



{-ir^^ (e^p-[^cothiK/M)yf''''^''y' 



F(M) dyf y - ^ - w ;..j ^^^ 

where we used the notations b± = Ke^''/*^/sinh (k/A/) and F(Af) = (M — 1)! for the Euler gamma- function. . 

Properties of this distribution arc discussed in much details in p6[ , also for the case of non-equivalent channels. 
Let us only briefly mention the most interesting features. 

Quick inspection of eq. (|l2) shows that it is indeed reduced to the x^ distribution, eq. (ph when the effective coupling 
to continua is weak: k <^ 1 . Under this condition resonances are typically too narrow to overlap with others: F ^ A. 
However, as long as the effective coupling becomes stronger: T — > 1, hence 5 S> 1 the parameter k grows large. Under 
these conditions another domain of resonance widths becomes more and more important: 

-K/Af sinh (k/M) K/Af sinh (k/M) 

"" {k/M) ^y^^^ (^/M) ' 

where the distribution eq.(p2[) shows the powerlaw decrease: p{ys) ~ 2^-^2/7^- The most drastic difference from 
eq.(1q) occurs for the maximal effective coupling g — I (i.e k — c»). In this regime the powerlaw tail extends up to 
infinity, making all positive moments ( starting from the first one) to be apparently divergent. One can argue that 
the powerlaw tail Myj"^ turns out to be dictated by classical processes of exponential escape typical for fully chaotic 
systems |16|. The rate of this escape in the semiclassical limit Af ^ 1 is provided by the value of the gap in the 
distribution of resonance width, see H for a more detailed discussion. 

The best candidates for checking the applicability of eq. ( |l^ ) to real physical systems are realistic models of ballistic 
mesoscopic devices subject to an applied magnetic field that serves to break the TRI pJ. It is however quite clear 
that all the basic qualitative features of the distribution eq.(|l3) (in particular, the powerlaw behavior p{ys) oc My~'^ 
for the overlapping resonance regime) should be valid for the systems with preserved TRI as well. Recent numerical 



data |57| support the vahdity of this conjecture. 

We have seen, that the non-Hermitian random matrix Hamitonian H^f appeared naturally in the scattering matrix 
description of open quantum systems. Actually, such a Hamiltonian is the most adequate tool to describe the quantum 
relaxation processes such as escape of the particle from the interaction region in the regime of quantum chaos. Some 
aspects of such a relaxation were studied some time ago in [ p8| and reconsidered in more details recently by Savin 
and Sokolov [M who used insights provided by the resonance widths distribution Eq.(n2h 

It is therefore quite natural that the resolvent of the non-Hermitian effective Hamiltonian H^f — H — iV is related 
to the probability for an excited system to decay via one of open channels. For this reason it can be used to calculate 
the statistics of such quantities as the Wigner time delay [ p9[ which is a measure of mean time spent by a scattered 
particle inside the interaction region. 

One more example of the utility of the resolvent of the non-Hermitian effective Hamiltonian Hef is that it can be 
related to the total photodissociation crossection in the regime of quantum chaos. The formula for the cross section 
autocorrelation function was derived recently by Fyodorov and Alhassid |23|] . Below we outline the derivation and 
present the final result. 

The total energy-dependent cross section a{E) is defined as a probability to be excited from a ground state \g) and 
to dissociate at a given energy E per unit time and per unit incoming photon fiux density. In the dipole approximation 
it is given by (see, e.g. detailed discussion in [gO[) : 

aiE)^aoJ2\{g\p\^T'{E))\' (13) 



Here jl is the dipole operator jl = e£r of the system in the external electric field £, ctq is a constant proportional 
to the excitation energy E and |$°"*(_B)) is the exact wave function of the system at energy E subject to outgoing 
boundary conditions in one of the open channels a = 1, 2, ..., M. 

It turns out, that for many systems of interest (e.g. molecules HO2 [^J] and H^ fljl), the cross section patterns 
look like irregular fluctuating signals consisting of many randomly positioned (partly) overlapping resonance peaks. 
Such a behavior ( typical also for the resonance scattering in atomic systems [ |l2y55|| ) has its origin in the underlying 
pattern of resonances in the complex plane. 

We already mentioned that one can relate the "internal" and "external" parts of the wavefunctions by using the 
resolvent of the non-Hermitian effective Hamiltonian H^f ■ In the present case such a relation can be writen as (see 
e.g. 0): 

E - Hlf) "' WB. (14) 

where M-component vector B contains amplitudes of outgoing waves in each of the open channels. 

The ground state wavefunction describes a bound state and as such has no components outside the interaction 
region. Using this fact and Eq.(|l4[) one finds after some algebraic manipulations that the cross section Eq.(|l3[) can 
be rewritten in the following convenient form: 



<j{E) ex Im ( g 



'^hu^' 



(15) 



which is just a particular case of the optical theorem. One also can arrive at the expression Eq.(n5|) by resumming 
the perturbation theory, see |B5[ for more details and relevant references. 

The advantage of the form Eq. (|5|) is that it expresses the photodissociation cross section in terms of the resolvent 
of an effective non-Hermitian operator 'He.f = Hin — ittWW'^ which is known to describe open chaotic systems in 
the random matrix formalism. It allows to apply very well developed methods of evaluating averages of products 
of resolvents based on the Efetov supermatrix formalism. Measuring energy separations in units of the mean level 
spacing of the closed system A one finds in such a calculation the cross section autocorelation function ||2^ : 

S{u = 27rf7/A) = j-^-^^ 1, 

to be a sum of two terms S(uj) — Si{uj) + 6*2 (w) given by the following expressions: 

S (u^^-CdxTdX rd\ ^°«KAiA2-A)](l~A^) f. (gg + A) 

(16) 



1 



where 



/i(A,Ai,A2)-(AiA2-A)2; hiX,Xi,X2) = 2X1X1 -Xf-Xl 



The parameters ga were introduced before and related to the transmission coefficients as ga — 2/Ta — 1. 

It is worth mentioning that each of the contributions S{lo)i^2 represent an interesting object by itself. Namely, 
Si{u>) coincides with the autocorrelation function of the so-called Wigner time delays studied in some detail in [p9[ , 
whereas 82(1^) is related by the Fourier-transform to the so-called "norm leakage" out the interaction region. The 
latter quantity was introduced recently by Savin and Sokolov as a characteristic of the process of quantum relaxation 
in chaotic systems and studied in detail for the simplest case of broken time-reversal invariance p4| . 

Actually, the starting Fermi golden rule formula Eq. ( [I4 ) is valid for an arbitrary excitation of the system with a 
weak perturbation fi. For this reason the autocorrelation function of crossections presented above is also of a general 
applicability. 

Finally, it is necessary to mention that in the limit Ta — for all a (corresponding to a closed system with purely 
bound spectra and no possibility for photodissociation) the expression Eq.(p^) reduces to the "oscillator strength" 
correlation function found by Taniguchi et al. ||62|. 



III. NON-HERMITIAN MATRICES WITH INDEPENDENT ELEMENTS: UNIVERSAL PROPERTIES 

IN THE REGIME OF WEAK NON-HERMITICITY 

To begin with, any N x N matrix J can be decomposed into a sum of its Hermitian and skew-Hermitian parts: 
J = Hi + iH2, where Hi — {J + J^)/2 and H2 = {J — J^)/2i. Following this, we consinder an ensemble of random 
N X N complex matrices J ^ Hi + ivH2 where Hp; p — 1,2 are both Hermitian: H^ = Hp. The parameter v is used 
to control the degree of non-Hermiticity. 

In turn, complex Hermitian matrices Hp can always be represented as Hi = Si + iuAi and H2 — S2 + iwA2, 
where Sp = Sp is a real symmetric matrix, and Ap = —Al is a real antisymmetric one. From this point of view the 
parameters u, w control the degree of being non-symmetric. 

Throughout the paper we consider the matrices Si, S2,Ai,A2 to be mutually statistically independent, with i.i.d. 
entries normalized in such a way that: 



lim —TrS?, = 


1 - -rr 

= lim —TrAr,At 


N^oo N P 


N^oo N P P 



1 (17) 

As is well-known [p|, this normalisation ensures that for any value of the parameter w ^ , such that u = 0(1) 
when A^ — > 00 statistics of real eigenvalues of the Hermitian matrix of the form H ^ S + iuA is identical (up to a 
trivial rescaling) to that of u = 1, the latter case known as the Gaussian Unitary Ensemble (GUE). On the other 
hand, for u = real eigenvalues of the real symmetric matrix S follow another pattern of the so-called Gaussian 
Orthogonal Ensemble (GOE). 

The non-trivial crossover between GUE and GOE types of statistical behavior happens on a scale u oc 1/N^'^ [ |63[ . 
This scaling can be easily understood by purely perturbative arguments [Q. Namely, for u c>c l/N^^^ a typical shift 
dX of eigenvalues of the symmetric matrix S due to the antisymmetric perturbation iuA is of the same order as the 
mean spacing A between unperturbed eigenvalues : JA ~ A ^ 1/A^. 

Similar perturbative arguments show pjq ], that the most interesting behavior of complex eigenvalues of non- 
Hermitian matrices should be expected for the parameter v being scaled in a similar way: v ex l/N^'"^ . It is just the 
regime when the imaginary part ImZk of a typical eigenvalue Zk due to non-Hermitian perturbation is of the same 
order as the mean spacing A between unperturbed real eigenvalues : ImZk ~ A ^^ 1/N. Under these conditions a 
non-Hermitian matrix J still "remembers" the statistics of its Hermitian part Hi. As will be clear afterwards, the 
parameter w should be kept of the order of unity in order to influence the statistics of the complex eigenvalues. 

It is just this regime of weak non-Hermiticity which we are interested in. Correspondingly, we scale the parameters 
asQ 

A. 

(18) 



2Vn' 2VN 

and consider a,(f),'w fixed of the order 0(1) when N —> 00. 

To be specific, we consider the real symmetric matrix Si to be taken from the ensemble of sparse random matrices 
[|65| , p6| characterized by the following probability density of a given entry Sij : 



^)S(S,,) + —. 



nSrj) = (1 - J:;)5{S,,) + j-MS^,) (19) 



where h{s) = h{—s) is an arbitrary even distribution function satisfying the conditions: /i(0) < c»; J h{s)s'^ds < 00 
and p stands for the mean value of non-zero matrix elements per column. Actually, this ensemble is the most general 
one among those with independent elements, and statistics of its eigenvalues was proved to be completely universal 
pa,p^, up to a rescaling by ensemble-dependent mean eigenvalue density i^{X). Statistics of the matrix elements of 
all other matrices 5*2,^1,2 is immaterial as long as their elements are statistically independent as well. 

The calculation of the mean density of complex eigenvalues follows essentially the same route as that outlined in 
the previous section for the resonances. The method used |M is a generalization of the Efetov's technique to the case 
of sparse random matrices suggested in ^^ (see some details also in [Q). As the result, one arrives at the following 
expression pO|: 



'^In the Letter |B9| there is a misprint in the definition of the parameter a. 



{p{X, y)) = ^^^ j d^i{Q)Str (ai^^q) Str (<7,g) exp ^S{Q) (20) 

2 9 l2 9 2 9 

5(g)=-^2/Str(<7,Q)-^Str(<7,g) + — Str (tsQ) - -^Str (^q) (21) 

where we introduced the scaled imaginary parts y = ttv{X)NY and used the notations: a? = {TTv{X)a) , 6^ ~ 
{■Kv{X)(j)) , (? = {Trh'(X)aw) . The supermatrices f2 and a entering this expressions are as foUows: 

T2 = diag{f3,f3}; '^ = [ f q 



and the supermatrices CTt and fs were defined after Eq.(|lO|). The expression (20) is just the universal a— model 
representation of the mean density of complex eigenvalues in the regime of weak non-Hermiticity we were looking for. 
The universality is clearly manifest: all the particular details about the ensembles entered only in the form of mean 
density of real eigenvalues i^{X). The density of complex eigenvalues turns out to be dependent on three parameters: 
a,b and c, controlling the degree of non-Hermiticity (a), and symmetry properties of the Hermitian part (6) and 
non-Hcrmitian part (c). 

The following comment is appropriate here. The derivation above can be done not only for ensembles with i.i.d. 
entries but also for any " rotationaly invariant" ensemble of real symmetric matrices Si . To do so one can employ 



the procedure invented by Hackenbroich and Weidenmiiller 1 68 1 allowing one to map the correlation functions of the 
invariant ensembles (plus perturbations) to that of Efetov's ct— model. 

Still, in order to get an explicit expression for the density of complex eigenvalues one has to evaluate the integral 
over the set of supermatrices Q. In general, it is an elaborate task due to complexity of that manifold. 

At the present moment such an evaluation was successfully performed for two important cases: those of almost- 
Hermitian matrices and real almost-symmetric 

matrices. The first case ( which is technically the simplest one) corresponds to (j) —^ c», that is 6 — > cxd. Under 
this condition only that part of the matrix Q which commutes with f2 provides a nonvanishing contribution. As the 

result, Str ( aQ ) = Str ( CTtQ ) so that second and fourth term in Eq.(pO|) can be combined together. Evaluating the 
resulting integral, and introducing the notation a? = a 



J--exp( — ^j / dtcosh{2ty)exp{-a'^t'^/2), 



Px(2/) = y-^exp(-^j / dtcosh{2ty)exp{-dH'/2), (22) 

where pxiv) is the density of the scaled imaginary parts y for those eigenvalues, whose real parts are situated around 
the point X of the spectrum (cf. Eq.(|lO|)). 

It is easy to see, that when a is large one can effectively put the upper boundary of integration in Eq.(E2p to be 
infinity due to the Gaussian cut-off of the integrand. This immediately results in the uniform density px{y) = (a^)~^ 
inside the interval \y\ < a^/2 and zero otherwise. Translating this result to the two-dimensional density of the original 
variables X,Y, we get: 

f . 2!t^ 2^ for irl < 27riy(X)v^(l + w^) 
p(X,Y) = i 47rt,2(i+t„2) I i_ V ; V ; /23) 

'^^ ^ \ otherwise ^ ^ 

This result is a natural generalization of the so-called "ellipic law" known for strongly non-Hcrmitian random 
matrices |^,^,^ . Indeed, the curve encircling the domain of the uniform eigenvalue density is an ellipse: ^^^^ , ^^v -t- 

1 as long as the mean eigenvalue density of the Hermitian counterpart is given by the semicircular law. The 
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semicircular density is known to be shared by ensembles with i.i.d. entries, provided the mean number p of non-zero 
elements per row grows with the matrix size as p ex N"; a > 0, see Q. In the general case of sparse or " rotationally 
invariant" ensembles the function ^{X) might be quite different from the semicircular law. Under these conditions 
Eq. (P3[) still provides us with the corresponding density of complex eigenvalues. 

The second nontrivial case for which the result is known explicitly is due to Efetov |32| . It is the limit of slightly 
asymmetric real matrices corresponding in the present notations to: cf) —f 0;w —^ oo in such a way that the product 
(j)w = c is kept fixed. The density of complex eigenvalues turns out to be given by: 

Px{y) = S{y) f diexp(-g2t2/2) + 2J-M H duexp f~^^C^^ f dttsmh{2t\y\)exp{-ch^ /2), (24) 



The first term in tiiis expression shows that everywhere in the regime of "weak asymmetry" c < oo a, finite fraction 
of eigenvalues remains on the real axis. 

Such a behavior is qualitatively different from that typical for the case of "weak non-Hermiticity" a < oo, where 
eigenvalues acquire a nonzero imaginary part with probability one. 

In the limit c >> 1 the portion of real eigenvalues behaves like c~^. Remembering the normalization of the 
parameter v, Eg. ([l7|) , it is easy to see that for the case of w = 0{1) the number of real eigenvalues should scale as 
\/N. The fact that of the order of N^^^ eigenvalues of strongly asymmetric real matrices stays real was first found 
numerically by Sommers et al. 1 29 il] , and proved by Edelman EM] . 



IV. GAUSSIAN ALMOST-HERMITIAN MATRICES: FROM WIGNER-DYSON TO GINIBRE 

EIGENVALUE STATISTICS. 



In the present section we concentrate on the particular case of almost-Hermitian random matrices with i.i.d. entries 
J = Hi +ivH2, where Hi and H2 are taken independently from the Gaussian Unitary Ensemble (GUE) of Hermitian 

matrices with the probability density ViX) = Q^^ exp [ - N/2J^ Tr xA , X ^ X"^ . 

Let us now introduce a new parameter r = (1 — -y^)/(l + u^) and choose the scale constant Jg to be equal to 
(1 + t)/2, for the sake of convenience. The parameter r controls the magnitude of correlation between Jj^. and Jkj- 
(JjkJkj) — t/X, hence the degree of non-Hermiticity. This is easily seen from the probability density function for our 
ensemble of the random matrices J: 



ViJ) = C],'exp 



N 



(l-r2) 



Tr(JJ^ - r Re P 



(25) 



where Cn = [tt (1 — t )/N ] ' . All the Jjk have zero mean and variance (| Jjfe| ) = 1/-/V and only Jjk and Jkj are 
pairwise correlated. If r = all the Jjk are mutually independent and we have maximum non-Hermiticity. When r 
approaches unity, Jjk and JL- are related via Jjk = J^,- and we are back to an ensemble of Hermitian matrices. 

Our first goal is to determine the n-eigenvalue correlation functions in the ensemble of random matrices specified 
by Eq. (Eq). The density of the joint distribution of eigenvalues in the ensemble is given by 



Vn{Zi 



.Z 



j^N{N+l)/2 



N 



r^l!---iV!(l-T2)^/2 



-AT ^ 



i=i 



r(^: 



z- 



}n 



z. 



Zkf. 



(26) 



j<k 



To derive Eq. (Eq) we integrate V{J) from Eq . (p5b over the surface of all complex matrices whose eigenvalues are 
Zi, . . . Zi\[. Following Dyson ||44[| , p. 501, see also |}48[| ) we decompose every complex matrix with distinct eigenvalues as 
J = U{Z + R)U'^ , where Z = diag{Zi, . . . Z^}, [/ is a unitary matrix, and i? is a strictly upper-triangular one. If we 
label the eigenvalues and require the first non-zero element in each column of U to be positive, then the decomposition 
is unique. The Jacobian of the transformation J —> {Z, R, U} depends only on Z and is given by the squared modulus 
of the Vandermonde determinant. So, integrating out R and U is straightforward and the resulting expression is 

Eq. (§)■ _ _ 

The form of the distribution Eq. ( Pq ) allows one to employ the powerful method of orthogonal polynomials [Q . 
Let Hn{z) denote the n-th Hermite polynomial. 



-«^^ 



cxpl^ 



dt t" exp ^izt 



The crucial observation borrowed from the paper J69| (see also the related paper [Q) is that the polynomials 




n = 0, 1, 2, . . ., are orthogonal in the complex plane Z = X + iY with the weight function 

N 



(27) 



(28) 



w (Z) = exp 



(l-r2) 



\Z\ 



-{Z^ + Z*^) 
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i. e. [ (P Zpni Z)pm( Z*)w'^ ( Z) = 5nm, where cPZ — dXdY. The standard machinery of the method of orthogonal 
polynomials ||44[| yields the n-eigenvalue correlation functions 



Rn{Zl, ..., Zn) — 

in the form 

Rn{Zl, . 

where the kernel Kpf(Zi, Z^) is given by 



Nl 



■ j d'Z^+l...(PZNVN{Z} 



{N -n)\ 
.,Zn)^dei[KN{Z„Zl)]l^^^,, 



(29) 



JV-l 



Kn{Zi,Z2) = 'w{Zi)w{Z^)^ Pn{Zi)pn(Z^ 



(30) 



With Eqs. (|2q)-(|30|) in hand, let us first examine the regime of strong non-Hermiticity, i.e. the case when 
limjv^oo(l — t) > 0. In this regime the averaged density of eigenvalues N^^Ri{Z) is asymptotically zero out- 
side the ellipse [Rez/(1 + t)]^ + [Im^/(1 - t)]^ < 1. Inside the ellipse \imN->oo N-^Ri{Z) = [tt{1 - r^)]"!. This 
sets a microscopic scale on which the averaged number of eigenvalues in any domain of unit area remains finite when 
N ^ oo. Remarkably, the r-dependence is essentially trivial on this scale: the statistical properties of eigenvalues are 
described by ^„(zi, . . . , 2:„) = iV-"i?„(\/]VZi, . . . , ^NZ,,) and 



lim i?„(zi,.. .,z„) 



1 



7r(l- 






det 



i,fc=i 



(31) 



This limiting relation can be inferred [ pO| from Mehler's formula for the Hermite polynomials W\. After the trivial 
additional rescaling z — > z\/\ — t^ the expression on the right-hand side in Eq. (|3l|) becomes identical to that found 
by Ginibre ^. 

Now we move on to the regime of weak non-Hcrmiticity. We know that in this regime new non-trivial correlations 
occur on the scale: ImZi.2 = 0{l/N), KeZi — R.CZ2 — 0{1/N). Correspondingly, we introduce new variables 
X, yi, 2/2, ^ in such a way that: x = Re {Zi + Z2)/2, yi.2 = iVIm (^1^2), ^ = -/VRe {Zi — Z2), and consider them finite 
when performing the limit iV — > cx). 

Substituting Eq.(p7|) into Eq.(|30[) and using the above definitions we can explicitly perform the limit N —^ 00, 
taking into account that liniAr^oo ^(1 ^ t) = a^/2. The details of the procedure are given elsewhere ^^. In this 
regime 



1 / uj/2 + iyi uj/2 + iy2 
lim — ttKm X H , X 



(32) 



1 



■ exp 



2/1 +2/2 , ix^yi -2/2) 



'^S.(X) ^^ 



TTUsc^x) 



exp 



"(2/1 + 2/2) + 2iuju 



with VsciX) = -^\/4 — X'^ standing for the Wigner semicircular density of real eigenvalues of the Hermitian part H 

of the matrices J. 

Equation ( p^ ) constitutes the most important result of the present section. The kernel K^ given by Eq. ( |3^ ) 
determines all the properties of complex eigenvalues in the regime of weak non-Hermiticity. For instance, the mean 
value of the density p{Z) = X]i=i '^^^^1^ ^ ^*) ^^ complex eigenvalues Z = X + iY is given by {p{Z)) = Kn{Z, Z*). 
Putting 2/1 = 2/2 and cj = in Eqs (|3j) we immediately recover the density Eq.(|22|) found by the supersymmetry 
approach 0. 

One of the most informative statistical measures of the spectral correlations is the 'connected' part of the two-point 
correlation function of eigenvalue densities: 



^In the present section we normalized H2 in such a way that for weak non-Hermiticity regime we liave limjv^oo Tr_ff| — N, 
whereas the nomalization Eq.([L7|) gives Iimjv_oo Tr_ff| — N{1 + w^). It is just because of this difference the parameter 5 
entering Eq.(p3) contains an extra factor 1 + w^ as compared to the present case. 
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{p{Z^)p{Z2))^ - {p{Z,)) 8'-^\Z^ - Z2) - 3^2(^1, ^2), 



(33) 



In particular, it determines the variance Y?{D) = {n{D)'^) — {n{D))'^ of the number n = Jjj(PZp{Z) of complex 
eigenvalues in any domain D in the complex plane as: 



MD) 



d'z, / d'Z2Mz^)piZ2)) - {piz^)){p{Z2))] 

D Jd 



<fZ{piZ)) 



d^Zi I d^Z2y2(Zi,Z2) (34) 



Comparing Eq.(p3|) with the definition Eqs. d2q)-(pfl) we see that the cluster function 3^2(^1,^2) is expressed in 
terms of the kernel Kn as 3^2(2'i, Z2) = \KN{Zi,Z*y^ 

It is evident that in the limit of weak non-Hermiticity the kernel Kj^ depends on X only via the semicircular density 
Vsc{X). Thus, it does not change with X on the local scale comparable with the mean spacing along the real axis 
A - 1/A^. 

The cluster function is given by the following explicit expression: 



3^(^,1/1,2/2) 



7V4 



du 



u(X) 



exp 



u{yi + 2/2) + iuu 



(35) 



The parameter a — Tiv{X)a controls the deviation from Hermiticity n When a — > the cluster function tends to 

GUE form 3^2(^,2/1, 2/2) = ^^(yi)(5(j/2) '^™ '^^■t — —■ In the opposite case a ^ 1 the limits of integration in Eq.(|35|) 
can be effectively put to ±00 due to the Gaussian cutoff of the integrand. The corresponding Gaussian integration 
is trivially performed yielding in the original variables Zi , Z2 the expression equivalent (up to a trivial rescaling) to 
that found by Ginibre ||: 3^2(^1,^2) = {N"^ /ira'^f eyi-p{-N'^\Zi - Z2\'^/a^}. 

The operation of calculating the Fourier transform of the cluster function over its arguments a;, j/1,2/2 amounts to 
simple Gaussian and exponential integrations. Performing them one finds the following expression for the spectral 
form-factor. 



b{qi,q2,k) 



duj 



dyi 



dy2y2{^, yi, 2/2) exp{27ri(cjfc + yiqi + 2/292)} 



(36) 



W4 / "' ^ 2 , 2 , o,.2M S'n[^"«'(gl+g2)('^(^)-|fc|)] ^, ,^, .... 

N exp{-— (gi +q2+2k )} ^ j-— ^ — - — '-9{i^{X) - \k\) 

2 Tf'a^{qi+ q2) 



where 0(u) = 1 for m > and zero otherwise. 

We see, that everywhere in the regime of weak non-Hermiticity < a < 00 the formfactor shows a kink-like 
behavior at |A:| = i^{X). This feature is inherited from the corresponding Hermitian counterpart-the Gaussian Unitary 
Ensemble. It reflects the oscillations of the cluster function with uj which is a manifestation of the long-ranged order 
in eigenvalue positions along the real axis |^ . When non-Hermiticity increases the oscillations become more and more 
damped. 

As we already discussed above the knowledge of the formfactor allows one to determine the variance S2 of a number 
of eigenvalues in any domain D of the complex plane. Small E2 is a signature of a tendency for levels to form a cristal- 
like structure with long correlations. In contrast, increase in the number variance signals about growing decorrelations 
of eigenvalues. 

In a general case this expression is not very transparent, however. For this reason we restrict ourselves to the 
simplest case, choosing the domain D to be the infinite strip of width Lx (in units of mean spacing along the real 
axis A = {i/sc{0)N)~^) oriented perpendicular to the real axis: < ReZ < L^A; —00 < ImZ < 00. Such a choice 
means that we look only at real parts of complex eigenvalues irrespective of their imaginary parts. It is motivated, 
in particular, by the reasons of comparison with the GUE case, for which the function T,{Lx) behaves at large L^ 
logarithmically: 'S{Lj.) oc InL^ ^. 

After simple calculations one finds n 



^2(Lx} — Lx 



'-^/ '§"-£' ='"^<'«« 



-(fr) 



(37) 



In our earlier Letter 
'in our earlier Letter fl 



B9| we used the definition of the parameter a difi'erent by a factor of 2 from the present one. 
M^the expression Eq.dSTf) and formulae derived from it erroneously contained tto instead of a. 
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First of all, it is evident that S2 grows systematically with increase in the degree of non-Hermiticity a = Tri'{Q)a. This 
fact signals on the gradual decorrelation of the real parts HeZi of complex eigenvalues. It can be easily understood 
because of increasing possibility for eigenvalues to avoid one another along the Y — ImZ direction, making their 
projections on the real axis X to be more independent. 

In order to study the difference from the Hermitian case in more detail let us consider again the large Lj: behavior. 
In that case the upper limit of the integral in Eq.(^7|) can be set to infinity. Then it is evident, that the number 
variance is only slightly modified by non-Hermiticity as long as a ^ Lx- We therefore consider the case a ^ 1 when 
we expect essential differences from the Hermitian case. 

In a large domain 1 <$^ L^ ^ a the second term in the integrand of Eg. (p7|) can be neglected and the number 
variance grows like 'E{Lx) — Lxf{Lx/a). We find it more transparent to rewrite the function f{u) in an equivalent 
form: 






dte 



which can be obtained from Eq.(p7[) after a simple transformation. 

For u = Lx/ a <^ 1 we have simply / ~ 1 and hence a linear growth of the number variance. For u ^ 1 we have 

/ « (^3/2u)-i. Thus, E2(L,) slows down: EaC^,) « (^. 

Only for exponentially large Lx such that In {Lx/ a) ^ a second term in Eq.(p7|) contributes significantly. Calculating 



its contribution explicitly and remembering that £3 |i 



ajijc'l'^) we finally find: 



^2{Lx > a) == 



-3/2 



A (... f %) - 1 



TT 



a 



where 7 is Euler's constant. This logarithmic growth of the number variance is reminiscent of that typical for real 
eigenvalues of the Hermitian matrices. 

Another important spectral characteristics which can be simply expressed in terms of the cluster function is the 
small-distance behavior of the nearest neighbor distance distribution |4jJ|,Q . 

We define the quantity p(Zo, S) as the probability density of the following event: i) There is exactly one eigenvalue 
at the point Z = Zq oi the complex plane, ii) Simultaneosly, there is exactly one eigenvalue on the circumference of 
the circle \Z — Zo\ = S iii) All other eigenvalues Zi are out of that circle: \Zi — Zo\ > S. 

As a consequence, the normalization condition is: / (PZq L dS p{Zq, S) = 1. In particular, for Hermitian matrices 
with real eigenvalues one has the relation: p{Zo,S) — S{lmZo)px{S), with pxiS) being the conventional "nearest 
neighbor spacing" distribution at the point X of the real axis Eil . 

We are interested in finding the leading small-S* behavior for the function p{Zo^ S). It turns out to be given by the 
following expression pO|: 



p{Zo,S) 



N 



de [{p{Zo)){p (Zo + ^e^'^)) - y2 {Zo, Zo + Se^')] 



(38) 



where we used the definition of the cluster function, Eq. ( p3| ) . 

In the regime of weak non-Hermiticity this formula is valid as long as the parameter S is small in comparison with 
a typical separation between real eigenvalues of the Hermitian counterpart: 5 <C A ~ 1/iV. 

Substituting the expression Eqs.(pl,p^ for the mean density and the cluster function into Eq.(|38|) one arrives after 
a simple algebra to the probability density to have one eigenvalue at the point Zq = X + iyoA and its closest neighbor 
at the distance \zi — zo\ = sA, A — {i'{X)N)~^, such that s ^ 1: 



pjX + iyoA,sA)\ 



s<.l 



TTiy^iX) 



'dyV 



yaiyo)-^:^ga{yo) 



-^—gaiya) 
oyo 



-i^S' 



^— / d9exp 
a^ Jo 



- — (s^cos^ 



2j/oscos( 



where 



9a{,y) 



du 



(27r)i/2 



exp{- 



2uy} 



First of all it is easy to see that in the limit a ^ 1 one has: Pa^i{Zo,s ^ 1) = -{s/a^)^ in agreement with the 
cubic repulsion generic for strongly non-Hermitian random matrices pSlEGlEfl]. On the other hand one can satisfy 
oneself that in the limit a — ^ we are back to the familiar GUE quadratic level repulsion: pa^o{Zo, s ^ 1) oc 6{yo)s'^. 
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In general, the expression Eq.(|3^) describes a smooth crossover between the two regimes, although for any a =/= the 
repulsion is always cubic for s — > 0. 

To this end, an interesting situation may occur when deviations from the Hermiticity are very weak: a <^ v2 and 
'observation points' Zq are situated sufficiently far from the real axis: 2|yo|/a ^ 2~^/^. 

Under this condition the following three regions for the parameter s should be distinguished: i)f ^ jf—^ iOirh' "^ 

± < 2l^ and finally iii) 2-^/2 < 2^^ <£<«-!. 



In the regimes i) and ii) the term linear in cos 6* in the exponent of Eq.(39) dominates yielding the result of 
integration to be the modified Bessel function tt/q (-fs^)- In the regime iii) the term quadratic in cos 9 dominates 

producing 27re~'^*/°^ /o [(s/a)^] « (27ra/s) . As the result, the distribution p{Zo, s) displays the following behavior: 

4- for - < 



PaiZa.s) 



d^ ( d 



9aiyo)T^goiyo) ~ ( -Kraoivo. 



Mva\ 

for T^ < 4 < 2%^ 



7r a 



for 2^ « ^ « a- 



with5o(2/) = 9a{v)\a=f)- 

Unfortunately, it might be a very difficult task to detect numerically the unusual power law p{s) ex s^' ^ because of 
the low density of complex eigenvalues in the observation points reflected by the presence of the Gaussian factor in 
the expression Eq. (p9|) . 

V. CONCLUSION 

In the present paper we addressed the issue of eigenvalue statistics of large weakly non-Hermitian matrices. 

Our original motivation came from the field of resonance chaotic scattering. The resonances, which are complex 
poles of the scattering matrix enter the theory as complex eigenvalues of a non-Hermitian effective Hamiltonian of 
a particular type: "Hef = H — iT. We demonstrated that one can extract mean density of such poles employing a 
mapping onto the supermatrix non-linear u— model. We also have shown how the resolvent of the non-Hermitian 
Hamiltonian Ti^f can be used to describe the process of chaotic photodissociation and presented the crossection 
autocorrelation function. 

Guided by our experience with the resonances, we found a regime of weak non- Hermiticity for other types of non- 
selfadjoint random matrices. The regime can be defined as that for which the imaginary part lu\Z of a typical complex 
eigenvalue is of the same order as the mean eigenvalue separation A for the corresponding Hermitian counterpart. 

Exploiting a mapping to the non-linear cr— model we are able to show that there are three different "pure" symmetry 
classes of weakly non-Hermitian matrices: i) almost Hermitian with complex entries ii) almost symmetric with real 
entries and iii) complex symmetric ones. Within each of these classes the eigenvalue statistics is universal in a sense 
that it is the same irrespective of the particular distribution of matrix entries up to an appropriate rescaling. There 
are also crossover regimes between all three classes. 

Our demonstration of universality was done explicitly for the density of complex eigenvalues of matrices with 
independent entries. Within the non-linear cr— model formalism one can easily provide a heuristic proof of such a 
universality for higher correlation functions as well as for "rotationally invariant" matrix ensembles, see p8[. The 
above feature is a great advantage of the supersymmetry technique. 

A weak point of that method is a very complicated representation of the ensuing quantities. It seems, that the 
explicit evaluation of the higher correlation functions is beyond our reach at the moment, and even a calculation of 
the mean density requires a lot of effort, see [ p8| , p2[ . As a result, at present time the mean density is known explicitly 
only for the cases i) and ii). 

Fortunately, because of the mentioned universality another strategy can be pursued. Namely, one can concentrate 
on the particular case of matrices with independent, Gaussian distributed entries for which alternative analytical 
techniques might be available. Such a strategy turned out to be a success for the simplest case of complex almost- 
Hermitian matrices, where we found the problem to be an exactly soluble one by the method of orthogonal polynomials. 



This fact allowed us to extract all the correlation functions in a mathematically rigorous way 1 39 5C[] . 

One might hope that combining the supersymmetric method and the method of orthogonal polynomials one will 
be able to elevate our understanding of properties of almost-Hermitian random matrices to the level typical for their 
Hermitian counterparts. 

From this point of view a detailed numerical investigation of different types of almost-Hermitian random matrices 
is highly desirable. Recently, an interesting work in this direction appeared motivated again by the theory of chaotic 



14 



scattering ||T^ . Unfortunately, matrices Tie/ emerging in that theory are different from the Gaussian matrices because 
of the specific form of the antihermitean perturbation iV necessary to ensure the unitarity of the scattering matrix. 
This fact makes impossible a quantitative comparison of our results with those obtained in pO| . The qualitative fact 
of increase in number variance with increase in non-Hcrmiticity agrees well with our findings. 
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